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Abstract 

Anomalous diffusion and non-Gaussian statistics are detected experimentally in a two- 
dimensional driven-dissipative system. A single-layer dusty plasma suspension with a Yukawa 
interaction and frictional dissipation is heated with laser radiation pressure to yield a structure 
with liquid ordering. Analyzing the time series for mean-square displacement (MSD), superdif- 
fusion is detected at a low but statistically significant level over a wide range of temperature. 
The probability distribution function (PDF) fits a Tsallis distribution, yielding q, a measure of 
non-extensivity for non-Gaussian statistics. 

PACS numbers: 52.27.Lw, 52.27.Gr, 82.70.Dd 
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Two-dimensional (2D) physical systems with a single-layer of mutually interacting parti- 
cles include a Wigner lattice of electrons on a liquid helium surface ions confined mag- 
netically in a Penning trap ^2], colloidal suspensions js], vortex arrays in the mixed state 
of type-II superconductors [4], and dusty plasmas levitated in a single layer [5,]. Within 
biological cell membranes, proteins can undergo 2D diffusion j^. 

We investigate random motion in driven-dissipative 2D liquids. In experiments with a 
si„g.e-layer dusty plasma suspension, particles interact with a Yukawa potential fl and are 
damped by gas friction. Like 2D granular media [8| , our particles are driven by an external 
source, and their average kinetic energy is set by a balance of driving forces and dissipation. 

Random motion is characterized using the mean-square displacement (MSD) time series 
and the probability distribution function (PDF). We define these separately for the x and 
y directions. For the x direction, MSD is {{x{At) — x{0))'^), and PDF is the histogram of 
|x(r) — x(0) I evaluated at a specified delay r. Random motion characterized as pure diffusion 



has two signatures: having a PDF that is Gaussian at all times 10| , and obeying the Einstein 
relation limf^oo((a;(At) — x(0))^) = 2D At [9] at long times, where D is a diffusion constant. 

Anomalous diffusion is a term used to describe motion that is found empirically to not 
exhibit these signatures. In this case. Pick's law is not satisfied, and there is no constant dif- 
fusion coefficient. Anomalous diffusion is associated with non-Gaussian statistics sometimes. 



as in a liquid with vortices 
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12l |. but not always [1 



Non-Gaussian statistics, which are often exhibited by non-equilibrium systems, can follow 
a Tsallis distribution for a random variable z, for example, displacement of a particle [l^ 



[l-(3{l-q)z 



211/(1-9) 



(1) 



where the non-extensivity parameter q is q = 1 for Gaussian and q 1 for non-Gaussian 
statistics. In addition to anomalous diffusion, Tsallis statistics has been used to describe 
turbulence .pure electron p,a..a.Q. 

There is a disagreement in the literature regarding whether motion can ever be purely 
diffusive in 2D liquids. There are simulations and theories predicting that motion is not 
diffusive for any type of interparticle interaction {igI . On the other hand, some recent 



simulations 



18 



19| have found signatures of diffusion in 2D liquids near the disordering 
transition, for various interparticle interactions. All these simulations can have limitations, 
including a finite system size. 
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The literature for diffusion in 2D liquids includes mostly simulations of simple liquids. 
Experiments with 2D liquids are fewer, including soft matter systems such as colloidal 



suspensions 
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22 



23 



24| . Distinguishing 



201 ]. granular material [2l|, and dusty plasmas 
diffusion from weakly anomalous diffusion at a statistically significant level requires large 
data sets for long times, which can be hard in a simulation and even harder in an 
experiment due to a limited thread lifetime. Here, the term "thread" refers to a time series 
of position measurements for one particle. 

Two-dimensional transport has been experimentally studied in dusty plasmas consisting 
of partially ionized rarefied gas that fills a 3D volume and contains charged microspheres 
suspended in a single thin layer. Juan et al. 23| measured the MSD transverse to a sheared 
flow in a quasi-2D multilayer suspension, and reported a diffusion coefficient that increased 
with the laser power that generated the shear. Nunomura et al. [Q] reported D for a 
layer that was perturbed by unstable movement of a few heavier particles in an incomplete 
second layer thereby raising the suspension kinetic temperature, T; they interpreted their 
results as indicating diffusion, with D cc T. Ratynskaia et al. 121] used a smaller cluster of 
microspheres suspended in a single layer that was heated presumably by instabilities in the 
dusty plasma; they described their motion as super diffusion, characterized by Levy statistics. 
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driven by flows 

shear la, ^51 and thermal conduction 



231 ]. There have also been other dusty plasma experiments to measure 
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271 ]. These experiments can be compared with 2D 



granular medium experiments, which are also driven-dissipative, and where diffusion takes 
a different character if there is a flow 28]. 

Here we report an experimental study of 2D diffusion with a different heating method. 
Our dusty plasma had no second layer {2^ . and we used a laser heating configured to avoid 
shear and macroscopic flows [l2, 23]. The method provides steady conditions, allowing us 



to record large data sets to provide a low detection limit for super diffusion. 

Our laser-heated suspension is a driven-dissipative system. Radiation pressure forces 
from rastered cw laser beams apply strong kicks in the ±x directions to a few particles at 
a time, and these collide with other particles to thermalize the motion. Energy input from 



the laser is ultimately balanced by dissipation on the neutral gas [29|, allowing a steady- 
state energy balance. Particles receive random kicks more often from collisions with other 
particles than from the laser beams. 

We used a radio-frequency (rf) apparatus [3^, Fig. [U^a). Plasma was formed from Ar at 
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8.6 mTorr by applying 13.56 MHz rf, yielding a dc self-bias of —95 V. We introduced > 6000 
polymer microspheres, 4.83 ± 0.08 fim in diameter, to form a single layer. The particles 
experienced a damping rate = 2.5 s^^ {2^. A cooled camera recorded 37 s movies of 
particle motion at 55.56 frames per second; its field-of-view (FOV) of 14.8 x 10.8 mm^ 
included about 900 particles in the central portion of the suspension. A 0.1 W Ar+ laser 
beam swept horizontally once per exposure illuminate particles. We verified that conditions 
remained steady throughout the experiment, and using a side-view camera we verified that 
no out-of-plane buckling occurred. 

We analyzed rnpvies, making > 2 x 10'' particle positions (x, y) using the moment method, 
optimized as in [311]. Tracking particles from frame to frame yields their trajectories, i.e., 
threads. Fig. [^b). Threads have a finite half-life T1/2, Fig. W^c), due to particles exiting 
the FOV or becoming indistinguishable during tracking. Lacking hydrodynamic coupling to 
any nearby wall, our suspension rotated as a rigid body when it was in a solid state. We 
always subtract this global rotation motion before calculating other quantities. We report 
temperatures and Ty computed from mean-square velocities. 



The Wigner-Seitz radius 32|] was a = 0.24 mm. The particle charge Q/e = —5700 and 
screening parameter a/A/) = 0.90 were measured using the natural phonon method 24 1. 
-Repeated measurements of Q and a/Xr) had a 4% dispersion. We calculated a time scale 



3^ ujj = 9.2 ms. 



To provide a well-controlled and adjustable particle kinetic temperature without affecting 
Q or a/Xo, we used the laser-heating method of 301]. A pair of 532 nm cw laser beams 
directed in the ±x directions drew Lissajous figures with frequencies fx = 40.451 Hz and 
fy = 25 Hz in a 16.2 mm wide rectangular stripe spanning the entire suspension in the x 
direction. We adjusted the temperature by varying the laser power. Fig. [^c). Gradually 
increasing the power of each beam to Pl = 0.75 W, as measured in the chamber, we melted 



the solid, as judged by Gg 
velocity. 



331 ] and g{r) 30|]. Flow velocities were < 8% of the thermal 



Compared to [30|], different Lissajous frequencies and a wider stripe helped reduce coher- 
ent modes and nonuniformities in the laser heating. Coherent modes, excited by the finite 
frequency of the laser-beam's repetitive motion, represented 10% of the total power for the 
X direction, and 6% for the y direction. We measure T and analyze random motion in each 
of six portions of the FOV, i.e., cells, in Fig. [T](b) inset; this method reduces the effect of a 
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35% spatial variation of temperature over the full FOV. Due to the anisotropic momentum 
input from our two laser beams, exceeded Ty, Fig. [T](c), with T^/Ty as large as 1.8. 

The MSD, Fig. [21 exhibits ballistic motion at short time, At < 0.1 s. For the x direction, 
the MSD also exhibits an artifact of laser heating between 0.1 and 1.0 s. Heating in the 
y direction arises only from collisions, so that it lacks this artifact, leading us to analyze 
motion mostly in the y direction. The MSD is always larger in the x direction, due to 
Tx > Ty. Here, we are interested in long-time behavior, characterized by measuring the 
diffusion exponent a [l^, which we present later. 

The PDF has a self-similar form, as demonstrated in Fig. [3] by the collapse of PDF 
curves at successive delays r ranging from 1 to 5 s. We scaled the displacement as = 
|y(r) — ?/(0)|/r°/^ to yield a self-similar form. 

A Tsallis distribution fits our self-similar PDF, and this fit yields q. Recall that Gaussian 
statistics would be indicated by g = 1. For the lower laser power in Fig. [3l we found 
q = 1.08 ± 0.01; for the higher power we found q = 1.05 ± 0.01. Repeating this analysis for 
other temperatures, and averaging the resulting q values within a temperature range, yields 
q in Table [11 

We observe q > 1, i.e., non-Gaussian statistics for all temperature ranges. These results 
are statistically significant, with a confidence level > 95% for rejecting the null hypothesis 
g < 1. The deviation from Gaussian statistics is greatest in our lowest temperature range, 
near the disordering transition. Unlike |l^, where a deviation from Gaussian statistics was 
attributed to prominent vortex flows, our experiment was prepared without flows. 

By calculating the diffusion exponent a, we determined that random motion exhibits 
superdiffusion. In this test, we fit MSD for displacements in the y direction to the scaling 
(At)", yielding a measurement of ay for different ranges of Ty and At. We then averaged the 
resulting ay values over ranges of Ty and time delay. In Table [H we report this mean value 
ay and a standard deviation of the mean. In all cases, we find ay > 1, i.e., superdiffusion, 
as compared with a = 1 reported in [2J]. The level of superdiffusion, ay — 1, is small, but 
statistically significant. Using our large data sets and the student's t-test, we can reject the 
null hypothesis that (Xy < 1 with a confidence level of > 95% (i.e., p < 0.05) in all cases 
except possibly the highest temperature range in Table [B 

The level of superdiffusion we detected did not reach the larger values of a > 1.1 observed 
in equilibrium 2D Yukawa simulations [l3|. Like our experiment, these simulations have no 
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TABLE I: The measure q of non-extensivity indicates non-Gaussian statistics. Mean diffusion 
exponent oty (and p-value for testing the null hypothesis that there is no superdiffusion) indicates 
super diffusion. 





time 




temperature range T,, (10'^ K) 






delay (s) 


10 - 20 


20-40 


40 - 60 


q 


1 < r < 5 


1.20 ±0.09 


1.11 ±0.06 


1.05 ±0.02 


y 


1 < At < 5 


1.052 ± 0.019 


1.059 ±0.011 


1.009 ± 0.011 


P 




0.007 


2 X 10"^ 


0.210 


ay 


5 < At < 9 


1.088 ± 0.048 


1.082 ±0.040 




p 




0.042 


0.038 




ay 


9 < At < 13 


1.146 ±0.062 






p 




0.028 






ay 


13 < At < 17 


1.183 ±0.064 






p 




0.006 







flows and no incomplete second layer, but they differ by assuming an energy equilibrium with 
no friction. We speculate that friction may suppress superdiffusion. Experimental systems 
experience friction, more strongly for colloids which are overdamped, and less strongly for 
dusty plasma, which uses a rarefled gas rather than a liquid solvent. 

Finally, we use our large data sets to investigate whether the MSD in our 2D liquid 
can yield a diffusion coefficient that has the scaling with temperature previously reported 
for an experiment 2^]. We tentatively calculate an estimated diffusion coefficient Dest as 
{{y{At) — ?/(0))^)/2At for displacement in the y direction, averaged for all movies and for 
various time delays At, yielding Fig. |H This coefficient of course increases with temperature. 



Based on a smaller data set, it was reported 



2J] that in a liquid this increase fits a straight 



line with an intercept on the horizontal axis at the melting temperature. However, we found 
that such a line flts our data more poorly {xl > 4) than a power law {xl = 0.98) or a 
straight line through the origin (x^ = 1.45). 

We thank Z. Donko, V. Nosenko, J. Scudder, and F. Skiff for helpful discussions. This 
work was supported by NASA and DOE. 
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FIG. 1: (a) Apparatus. An argon plasma is formed in a vacuum chamber, not shown here. Micro- 
spheres suspended in a single layer above a lower electrode are kicked by a pair of rastered laser 
beams, (b) Particle trajectories for = 2.5 W in cell 1. We divide the camera's full field-of-view 
(FOV) into 6 cells, as shown in the inset. Motion is mostly ballistic in the 0.5 s time interval shown 
here, (c) Heating. Kinetic temperatures, averaged over the full FOV, increase with P^. Structures 
with Ge < 0.1 were identified as liquid and lacked large crystalline domains; structures identified 
as solid were not defect-free. At higher temperatures, threads have a shorter half- life tx/2- The 
anisotropy > Ty is due to laser heating. 



FIG. 2: (Color online) MSD, averaged over threads and overlapping time segments [3J], computed 
separately for x and y. Data here are for two laser powers. The MSD in the x direction is larger 
than that in the y direction, due to Tx > Ty. Distance is normalized by the Wigner-Seitz radius a. 



FIG. 3: (Color online) PDF for the y direction, and the same conditions as in Fig. [2j Self- 
similarity is revealed by the collapse of PDFs at successive delays r, with the scaled displacement 
AC = |y(r) - y(0)|/r°/2. Fitting to Eq. ^ yields the smooth solid curves and the measure of 
non-extensivity q. Other curves shown bracket the smooth curve for the fit. Data here are for 
1 < r < 5 s. 



FIG. 4: (Color online) Temperature dependence of D^st calculated from {{y{At) — ?/(0))^)/2At. 
The smooth curve is a fit, weighted by error bars, of data for a liquid. Scatter is due to finite time 
interval and number of threads. 
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This figure "Figure3_sep_ll_2007_4.gif" is available in "gif" format from: 
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